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Abstract 

Quasi Elastic Scattering (QEL) is the dominant mechanism for 
producing both Fully Contained Events and Partially Contained Events 
in the SuperKamiokande (SK) detector for atmospheric neutrinos in 
the range ~0.1 GeV to '^lO GeV. In the analysis of SK events, it is 
assumed that the zenith angle of the incident neutrino is the same as 
that of the detected charged lepton. 

In the present paper, we derive the distribution function for the 
scattering angle of charged leptons, the average scattering angles and 
their standard deviation due to QEL. Thus, it is shown that the SK 
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assumption for the scattering angle of the charged lepton in QEL is not 
valid. Further, for the rigorous analysis of the experimental data, the 
correspondence between the zenith angle of the charged leptons and 
that of the incident neutrinos should be examined by taking account 
of the influence of the azimuthal angle of the charged particle on its 
zenith angle. We conclude that it is not possible to reliably assign the 
zenith angle of the incident neutrino to that of the charged lepton, 
and so the zenith angle distribution of charged leptons under the SK 
assumption does not reflect that of the real incident neutrinos. 

This result has clear implication for attempts to detect neutrino 
oscillations from the analyses of Fully Contained Events and Partially 
Contained Events in SuperKamiokande. 

1 Introduction 

The report of oscillations between muon and tau neutrinos for atmospheric neu- 
trinos detected with SuperKamiokande (SK, hereafter) is claimed to be robustly 
established for the following reasons: 

(1) The discrimination between electrons and muons in the SK energy range, 
say, ~0.1 GeV to ~1 GeV, has been proved to be almost perfect, as demonstrated 
by calibration using accelerator beams p. 

(2) The analysis, for both Fully Contained Events and Partially Contained 
Events, of the zenith angle distribution of electron-like (single ring) events and 
muon-like (single ring) events, based on the well established discrimination pro- 
cedure mentioned in (1), reveals a significant deficit of muon-like events but the 
expected level of electron-like events. It is concluded that muon neutrinos oscil- 
late into tau neutrinos which cannot be detected due to the small geometry of SK 
[21. The analysis of SK data presently yields sin^26' > 0.89 and 1.8 x lO^^eV^ < 
Am^ < 4.5 X 10"^eV^ at 90% confidence level. ^ 

(3) The analysis of Upward Through Coing Events and Stopping Events, in 
which the neutrino interactions occur outside the SK detector, leads to similar 
results to (2). The charged leptons which are produced in these categories are 
regarded as being exclusively muons, because electrons have negligible probabilities 
to produce such events as they lose energy very rapidly in the sorrounding rock. 

Thus for these events the discrimination procedure described in (1) is not 
required, and , therefore, the analysis here is independent of the analysis in (2). 



^The numerical values of the combination of sin 29 and 5m^ change from time to time. 
Here, we cite those from the most recent work j2j. 



For these events, however, the SK group obtains the same parameters for neutrino 
oscihations as in (2). ^ 

The most robust evidence for neutrino oscillations is regarded as being that 
from the analysis of Fully Contained Events, because all details of the events 
are measured inside the detector and its analysis is free of the ambiguities that 
arise in the other analyses. Partially Contained Events have uncertainties due to 
their unknown ending points, and Stopping Muon Events and Upward Through- 
Coing Muon Events both have uncertainties due to the unknown starting points, 
preventing a complete analysis of the events concerned. 

In the SK experiment, the direction of the incident neutrino is assumed to 
be the same as that of the emitted charged lepton, i.e., the (anti-)muon or (anti- 
)electron, in the analysis of both Fully Contained Events and Partially Contained 
Events |^. However, as we will show, this assumption does not hold in the 
most important energy region, ~0.1 GeV to ~10 GeV, for both Fully Contained 
Events and Partially Contained Events, where Quasi Elastic Scattering (QEL) of 
the neutrino interaction is most dominant 0. In the present paper, we examine 
the implication of this assumption in a quantitative way. 

2 Cross Sections of Quasi Elastic Scattering 
in the Neutrino Reaction and the Scatter- 
ing Angle of Charged Leptons. 

We examine the following reactions due to the charged current interaction (c.c). 

f e + n — > p + e~ 
v^ + n — > p + fJ.~ 

^e+ P — ^ n + e"*" (1) 

^^l+P — > n + /i"*" 

The differential cross section for QEL is given as follows [S]. 



da _G>s%/^^^2^^5(q2) z_^ ^c{Q'^ z__L I (2) 



dQ2 8ttE^ 



's — u 

[ M2 J 


+ C{Q^) 


' s — u 
. M2 



^It seems strange for different quality experimental data to give similarly precise results, 
because Fully Contained Events are of higher experimental quality compared with those 
of both Upward Through Going Muon and Stopping Muon Events. See also footnote 1. 



where 



A{Q^ 
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^(£-A^f.f'X^f^(^-^)^a^(4^^) 



^'[w-V^^'^'Tp^^'[w-4mr^'v^My 



B{Q') = {fi + f2)9iQ' 



4 V ' '4Af2 

The signs + and — refer to I'fj.fe) and i^fj,(e) for charged current (c.c.) interactions, 
respectively. The Q^ denotes the four momentum transfer between the incident 
neutrino and the charged lepton. Details of other symbols are given in jS]. 

The relation among Q^ , E^, the energy of the incident neutrino, Ei, the energy 
of the emitted charged lepton (muon or electron or their anti-particles) and ^s, the 
scattering angle of the charged lepton, is given as 

Q^ = 2E^E(,{l-coses). (3) 

Also, the energy of the charged lepton is given by 

E, = K- ^. (4) 

In Figs. 1 to 4, we give the differential cross sections for different charged 
leptons and for different incident neutrino energies. It is clear from these figures 
that the cross sections of the neutrinos are larger than those of anti-neutrinos in 
the lower energy region, say, ~0.1 GeV to ~10 GeV and that their difference are 
negligible in higher energy region, say, beyond ~10 GeV. The differences between 
(anti-)muon and (anti-) electron are negligible except in the lower energy region, 
below ~0.1 GeV. 

In Fig. 5, we give the total cross sections for muon neutrinos and anti- muon 
neutrinos. It can be seen that the differences between thein are rather large in the 
lower energy region, say, ~0.1 GeV to ~10 GeV, but that the cross sections have 
similar values above ~10 GeV. The corresponding total cross sections for electron 
and positron are almost the same as those of the muon neutrinos and anti-muon 
neutrinos, respectively. 

Now, let us examine the magnitude of the scattering angle of the charged lepton 
in a quantitative way, as this plays a decisive role in determining the accuracy of the 
direction of the incident neutrino, which is directly related to the reliability of the 
zenith angle distribution of both Fully Contained Events and Partially Contained 
Events in SK. 
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Figure 1: Differential cross section for QEL for muon neutrinos witli different 
incident energies, 0.5 GeV, IGeV and 5 GeV. 
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Figure 2: Differential cross section for QEL for anti-muon neutrinos with 
the same incident energies as in Fig. 1. 
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Figure 3: Differential cross section for QEL for electron neutrinos with 
diffrent incident energies, 0.1 GeV, 1 GeV, 10 GeV and 100 GeV. 
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Figure 4: Differential cross section for QEL for anti-electron neutrinos with 
the same incident energies as in Fig. 3. 
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It can be seen from Eqs. (3) and (4) that there is a relation between the 
energy Ei of the charged lepton and its scattering angle Og for a given incident 
neutrino energy E^. Figure 6 shows this relation for muon, from which we can 
easily understand that the scattering angle ^s of the charged lepton ( muon here 
) cannot be neglected. For a quantitative examination of the scattering angle, we 
construct the distribution function for 9s of the charged lepton from Eqs. (2) to 
(4) by using the Monte Carlo method. '^ 

Figure 7 gives the distribution function for ^s of the muon produced in the 
muon neutrino interaction. It can be seen that the muons produced from lower 
energy neutrinos are scattered over wider angles and that a considerable part of 
them are scattered even in backward directions. Similar results are obtained for 
anti-muon neutrinos, electron neutrinos and anti-electron neutrinos. 

Also, in a similar manner, we obtain not only the distribution function for the 
scattering angle of the charged leptons, but also their average values < 9s > and 
their standard deviations 0"=. Table 1 shows them for muon neutrinos, anti-muon 



''The distribution functions for the scattering angle can be more easily obtained by the 
numerical method than done by the Monte Carlo method. The reason why we obtain them 
by the Monte Carlo method is to apply this Monte Carlo simulation procedure to another 
purpose by which we get the essential part of this paper. The distribution function for 
scattering angle in the inelastic scattering of the neutrino interaction and other quantities 
had been obtained by Kobayakawa et al. See |7]. 
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Figure 6: Relation between the energy of the muon and its scattering angle 
for different incident muon neutrino energies, 0.5 GeV, 1 GeV, 2 GeV, 5 GeV, 
10 GeV and 100 GeV. 
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Figure 7: Distribution functions for the scattering angle of the muon for 
muon-neutrino with incident energies, 0.5 GeV, 1.0 GeV and 2 GeV. Each 
curve is obtained by the Monte Carlo method (one million sampling per each 
curve). 



Table 1: The average values < ^s > for scattering angle of the emitted 
charged leptons and their standard deviations as for various primary neutrino 
energies E^. 
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neutrinos, electron neutrinos and anti-electron neutrinos. In the SK analysis, it is 
assumed that the scattering angle of the charged particle is zero (31 15 > hut Table 1 
demonstrates that such an assumption does not hold at all in both Fully Contained 
Events and Partially Contained Events. 

Therefore, it is surely concluded that the zenith angle distribution obtained by 
SK for Fully Contained Events and Partially Contained Events are almost unreli- 
able, which will be shown later more concretely. 



3 Influence of Azimuthal Angle of Quasi Elas- 
tic Scattering over the Zenith Angle of the 
Fully Contained Events 

For three incident cases(vertical, horizontal and diagonal), Fig. 8 gives a schematic 
representation of the relationship among the zenith angle of the incident neutrino, 
Oy, the scattering angle of the charged lepton, ^s; and the azimuthal angle of the 
charged lepton, </>. In this paper, we measure 6y from the vertical upward direction. 
From Fig. 8-a, it can been seen that the zenith angle 9 of the charged lepton 
is not influenced by its (j) in the vertical incidence {Ou = 0°). From Fig. 8-b, 
however, it is obvious that the angle influences on its zenith angle greatly in 
the case of horizontal incidence of the neutrino {9^ = 90°). Namely, half of the 
charged leptons are recognized as upward going, while the other half is classified as 
downward going. The intermediate case (diagonal incidence of 9,^ = 43°) between 
the above two cases is shown in Fig. 8-c. As discussed above, we must take 
account of the effect due to the azimuthal angle (j) in the analysis of both Fully 
Contained Events and Partially Contained Events. Figure 9 shows the relation 
between direction cosines of the incident neutrinos, {£,m,n), and those of the 
charged lepton, (4,mr,nr), for certain 9s and <p. Then, we have the following 
expression: 
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' sin^scosc; 
sin^sSin(? 
\ cos6's. 



(5) 



where n = cos9y, and n^ = cos9. Here, 9 is the zenith angle of the charged lepton. 
Using Eq. (5), we carry out a Monte Carlo calculation to examine the influence 
oi (j) on 9 for the charged lepton. The Monte Carlo procedure for the determination 
of the real 9 of the charged lepton whose parent (anti-)neutrino has fixed 9i, and 
E^ involves the following steps: 

1. We extract Q^ from the probability function for the differential cross section 
with a given Ei, (Eq. (2)) by the random sampling. 

2. We obtain Ei from Eq. (4). 

3. We obtain 9s from Eq. (3). 

4. We decide d>, which is obtained from 



2-K^. 



(6) 
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and those of the emitted charged lepton. 



12 



Here, .^ is a uniform random number between and 1. 

5. We obtain {£r,mr,nr) by using Eq. (5) . The n^. is the cosine of the zenith 
angle of the charged lepton which should be contrasted to n, that of the incident 
neutrino. In the SK analysis, it is assumed that 

{ir,mj.,nr) = {i,m,n). (7) 



We examine 6, the zenith angles of the charged muons, for three typical zenith 
angles of the incident neutrino with Ei^ = 1 GeV, namely, cos 9^ = 1 (vertically 
upward), cos9,y = 0.731, i.e., 9,/ = 43° (diagonal in SK), and cosO^ = (horizon- 
tal). 

A : vertically upward incident neutrino events 
In this case, as easily understood from Fig. 8-a, we can skip the step 4 mentioned 
above, because the change of (j) of the charged muon has no influence over its zenith 
angle 9. 

Here, although the direction of the charged muon approximately retains the 
primary direction of the incident neutrino at higher energies, say, above ~20 GeV, 
it deviates from the primary direction (vertical), even to backward directions, at 
lower energies, say, below ~1 GeV (see also Fig. 7). 

In Fig. 10, we show the scatter plot between E^, the emitted energy of the 
charged muon and cos 6, the cosine of the zenith angle of the charged muon for 
the vertically incident neutrino. In this case, the relation between 9 and Ei, is 
unique because of independence of the azimuthal angle. The reason is as follows: 
For a given Q^ (see Eq.(2)) we obtain the energy of the charged lepton and its 
scattering angle uniquely due to the two body kinematics and the zenith angle of 
the charged lepton are never influenced from their azimuthal angle, because 9 is 
measured from the vertical direction , parallel to the axis of the detector. 

B: horizontally incident neutrino events 
This case shows a great contrast to the vertical case. In Fig. 11, we give the scatter 
plot between E^y and cos^.In this case, the influence of (p over 9 is the strongest 
amomg the three, because this influence appears through Eq.(5) as strongly as pos- 
sible. The cos9 has the widest distribution, even if the charged muon has the same 
energy. As 9 is symmetrical due to its azimuthal symmetry, the scattered muons 
have almost lost the direction of primary neutrinos at lower energies, say, ~1 GeV 
and show a uniform distribution between upward and downward. It should be 
noted that the scatter plots are distributed symmetrically around cos 9 = 0, re- 
flecting the azimuthal symmetry. 
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Figure 10: The scatter plot between the fractional energies of the produced 
muons and their zenith angles for vertically incident muon neutrinos with 
1 GeV. The sampling number is 1000. 
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Figure 11: The scatter plot between the fractional energies of the produced 
muons and their zenith angles for horizontally incident muon neutrinos with 
1 GeV. The sampling number is 1000. 
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14 



C: diagonally upward incident neutrino events 
In Fig. 12, we show the similar scatter plot for the diagonal case, which shows the 
intermediate situation between Figs. 10 and 11, as expected. From Fig. 12 for 
01, = 43°, it is still apparent that (j) influences 6* in a considerable degree, even if 
the energies of the charged muon are the same. 

We can get the similar scatter plots for electrons to those for muons. 

We sum the events concerned over E^ for a given cos9 in Figs 10 to 12 and 
show the results of the zenith angle disributin in Figs 13 to 15. In Fig. 13, the 
neutrino enters vertically. In this case, even if we define cos 9 = 0.8 ~ 1.0 as 
"vertical" , about one half of the real events are not recognized as vertical. In Fig. 
14 we give the zenith angle distribution for the horizontal(cos ^,^ = 0). Comparing 
it with Fig. 13, it is easily understood that the direction of the primary neutrino 
here has been almost completely lost, namely, nr is very different from n. Thus, 
the SK assumption {rij- = n) clearly does not hold in this energy region. In Fig. 
15, we show the intermediate state between Figs. 13 and 14. The tendencies in 
Figs. 13 to 15 for muon neutrinos are quite similar to those for electron neutrinos. 

Furthermore, the dependence of cos9 on cos^jy becomes more sensitive at lower 
energies, say, below ~0.2 GeV. Also, they are a little weaker for anti-neutrinos 
than for neutrinos, as can be inferred from the cross sections in Figs. 1 to 4 and 
Fig. 5. 

It should be, here, emphasized that the azimuthal angle of the charged lepton 
greatly influences discrimination between Fully Contained Events and Partially 
Contained Events, which is strongly dependent on the generation points of the 
neutrino events concerned. The estimation of this discrimination is very impor- 
tant for getting Fully Contained Events, because the correct analysis of Fully Con- 
tained Events only in the SK experiment makes it possible to lead less ambigious 
conclusion on neutrino oscillation. @ 

4 The Zenith Angle of Fully Contained Events 
Taking Account of the Overall Neutrino 
Spectrum 

In the previous section, we obtain the zenith angle distribution of both Fully Con- 
tained Events and Partially Contained Events for a given energy of the incident 
neutrino, taking account of the effect of the azimuthal angle for a given zenith 
angle of the incident neutrino. However, in order to examine the effect of 9s and 
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(j) of the charged lepton in QEL on the direction of primary neutrinos, we must 
consider the overah neutrino spectrum, not a fixed neutrino energy. 

The concrete procedures for a selected event with a given O^, of the incident 
neutrino are summarized again as follows: 

Procedure A 
We formulate the differential neutrino interaction probability functions for muon 
(electron) and anti-muon (anti-electron) in which the (anti-)neutrino differential 
energy spectra are combined with the total cross sections of the corresponding 
neutrino concerned for QEL. In such a calculation, the differential neutrino spectra 
at Kamioka site which cover from 0.1 GeV to 100 GeV obtained by Fiorentini et 
a/.jHl are utilized. 

Procedure B 
The determination as to whether the (anti-)neutrino concerned belongs to either 
neutrino or to anti-neutrino is done by random sampling according to the proba- 
bility functions constructed from Procedure A. 

Procedure C 
We randomly sample the energy of the (anti-)neutrino concerned according to the 
probability function constructed from the Procedure A. 

Procedure D 
We randomly sample Q^ from Eq. (2) whose neutrino energy is determined by 
Procedure C. 

Procedure E 
We decide the energy of the charged lepton concerned, E^, for the Q^ from Eq. (4). 

Procedure F 
We decide the scattering angle of the charged lepton concerned, 9^-, from Eq. (3). 

Procedure G 
We randomly sample the azimuthal angle of the charged lepton concerned, 0, from 
Eq. (6). 

Procedure H 
We decide the direction cosines of the charged lepton concerned by using Eq. (5), 
being accompanied by Procedures F and G. 
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We repeat Procedures B to H until we reach the desired trial number. 

In Figs. 16 to 18, we give the zenith angle distributions of the sum of /i+ and 
/i~ for a given zenith angle of parent neutrinos. In Figs. 19 to 21, we give those 
for 6+ and e~ which correspond to Figs. 16 to 18. There are no major differences 
between (anti-)muon and (anti-)electron. The small differences between them are 
mainly due to the different neutrino fluxes between (anti-)muon neutrinos and 
(anti-)electron neutrinos. 

If the SK assumption is valid, then, the zenith angle distribution of the charged 
leptons should be of the delta function type. However, the real distributions of 
the charged leptons in Figs. 16 to 21 are quite far from the delta function type 
distribution. From Figs. 16 to 21, we can conclude that the present SK analysis 
procedures for Fully Contained Events and Partially Contained Events can not 
essentially determine the zenith angles of the incident neutrinos. It should be 
finally noticed that we must consider not only the contribution of the downward 
neutrino events which originate from the upward going neutrinos, but also that of 
the upward neutrino events which originate from the downward going neutrinos. 

5 Conclusions 

The standard SK analysis for Fully Contained Events and Partially Contained 
Events assumes that the direction of the emitted charged lepton is the same as 
that of the incident neutrino. However, Figs. 16 to 22 clearly show that the 
zenith angle distributions of the charged leptons have a wide spread for a fixed 
zenith angle of the incident (anti-)neutrino, with the leptons even being scattered 
backward in some cases, while the corresponding distributions obtained from the 
SK assumption should be of the delta function type, the peak of which correspond 
to the zenith angles of the incident neutrinos. 

Thus, it is clearly shown that the real zenith angle distributiuon of the charged 
particle in which their scattering angle is correctly taken into account cannot be 
approximated in any sense to the delta function type distribution which comes 
directly from the SK assumption. 

Consequently, the SK assumption (n^. = n) does not reflect the real zenith 
angle for directions of the incident neutrinos. 

Namely, our results show clearly that the zenith angle distribution of the muon- 
like (single ring) events and electron-like (single ring) events cannot reflect that of 
incident neutrinos, even if the discrimination between electron and muon is perfect. 
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Figure 16: Zenith angle distribution of /i^ and fi^ for z/^ and z/^ with the ver- 
tical direction, taking account of the overall neutrino spectrum at Kamioka 
site. The sampling number is 10000. 
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Figure 17: Zenith angle distribution of /i~ and /i"*" for z/^ and u^ with the hor- 
izontal direction, taking account of the overall neutrino spectrum at Kamioka 
site. The sampling number is 10000. 
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Figure 18: Zenith angle distribution of fi~ and //+ for z/^ and z/^ with the di- 
agonal direction, taking account of the overall neutrino spectrum at Kamioka 
site. The sampling number is 10000. 
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Figure 19: Zenith angle distribution of e~ and e"'" for i^e and z/g with the ver- 
tical direction, taking account of the overall neutrino spectrum at Kamioka 
site. The sampling number is 10000. 
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Figure 20: Zenith angle distribution of e~ and e^ for z/g and z/g with the 
horizontal, taking account of the overall neutrino spectrum at Kamioka site. 
The sampling number is 10000. 
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Figure 21: Zenith angle distribution of e^ and e+ for z/g and z/g with the 
diagonal, taking account of the overall neutrino spectrum at Kamioka site. 
The sampling number is 10000. 
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The zenith angle distribution of both Fully Contained Events and Partially 
Contained Events in our numerical computer experiment in which the polar angle 
(the scattering angle) as well as the azimuthal angle of the charged leptons are 
correctly taken into account will be published elsewhere. 



^We have studied the SK discrimination between electrons and muons by simulating 
electron-like (single ring) events and muon-like (single ring) events and examining the 
discrimination between them, based on the concept of pattern recognition. We conclude 
that the level of discrimination claimed by SK is practically impossible to achieve. This 
study will be described in a later publication. 
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